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We develop perturbative QCD formalism for three-body nonleptonic B meson decays. Leading 
contributions are identified by defining power counting rules for various topologies of amplitudes. 
The analysis is simplified into the one for two-body decays by introducing two-meson distribution 
amplitudes. This formalism predicts both nonresonant and resonant contributions, and can be 
generalized to baryonic decays. 

rn"; 

' The fundamental concept of perturbative QCD (PQCD) is to separate hard and soft dynamics in a QCD process. 

The former is calculable in perturbation theory, while the latter, though not calculable, is treated as a universal input. 
' The separation can be performed in the framework of collinear factorization [1] or of Ut factorization [2,3], in which 
an amplitude is expressed as a convolution of a hard kernel H with a hardon distribution amplitude $(a;) or with 
a hadron wave function $(a;,fcT)j x and kx being a longitudinal momentum fraction and a transverse momentum, 
respectively. Collinear factorization works, if it does not develop an end-point singularity from a; — > in the above 
convolution. If it does, collinear factorization breaks down, and fcy factorization is more appropriate. 
' It has been known that collinear factorization of charmed and charmless two-body B meson decays suffers end- 
point singularities. The PQCD formalism for these modes based on kx factorization theorem was then derived [4-6], 
which has been shown to be infrared-finite, gauge- invariant, and consistent with the factorization assumption in 
^ , the heavy-quark limit [7,8]. If one still employs collinear factorization, an alternative approach, the so-called QCD- 
£2 ' improved factorization [9], can be developed. In this approach the end-point singularities in the leading contributions 
, are absorbed into B meson transition form factors, and those appearing at the subleading level are regularized by 
■ arbitrary (nonuniversal) infrared cutoffs of momentum fractions x. Without the arbitrary cutoffs, PQCD has a 
\ predictive power, whose predictions for B PP, VP, and VV modes are all in agreement with data [10]. 
fS| ' Three-body nonleptonic B meson decays have been observed recently [11,12]. Viewing the experimental progress, it 
' is urgent to construct a corresponding framework. Motivated by its theoretical self-consistency and phenomenological 
(— I , success, we shall generalize PQCD to these modes. A direct evaluation of the hard kernels, which contain two virtual 
O i' gluons at lowest order, is not practical due to the enormous number of diagrams. On the other hand, the region 
' , with the two gluons being hard simultaneously is power-suppressed and not important. Therefore, a new input is 
necessary in order to catch dominant contributions to three-body decays in a simple manner. The idea is to introduce 
(-H two-meson distribution amplitudes [13], by means of which a factorization formula for a _B — s- /ii/i2^3 decay amplitude 
• • is written, in general, as 

! M = <i>B^ H (E) <^h^h2 ® */i3 • (1) 

^ ' 

It will be shown that both nonresonant contributions and resonant contributions through two-body channels can be 
included through the parametrization of the two- meson distribution amplitude $^1/12- 

Three-body decay amplitudes are classified into four topologies, depending on number of light mesons emitted from 
the four-fermion vertices. Topologies I and III, shown in Figs. 1(a) and 1(c), are associated with one light meson 
emission and three light meson emission, respectively. The bubbles denote the distribution amplitudes, which absorb 
nonpcrturbative dynamics. The hard kernel H contains only a single hard gluon exchange. The former involves 
transition of the B meson into two light mesons. In the latter case a B meson annihilates completely. For two light 
meson emission shown in Fig. 1(b), we assign lis to the special amplitude corresponding to the scalar vertex, and II 
to the rest of the amplitudes. Both topologies II and lis are expressed as a product of a heavy-to-light form factor 
and a time- like light- light form factor in the heavy-quark limit. 

The dominant kinematic region for three-body B meson decays is the one, where at least one pair of light mesons has 
the invariant mass of 0{AMb) for nonresonant contributions and of O(A^) for resonant contributions, A — Mb — nib 
being the B meson and b quark mass difference. An example is the configuration, where all three mesons carry 
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momenta of 0{Mb), but two of them move almost parallcUy. In the above dominant region coUinear factorization 
theorem apphes to topology I, since it is free of end-point singularities as shown below. With the pair of mesons 
emitted with a small invariant mass, the evaluation of topologies II and lis is the same as of two-body decays. 
The contribution from the region, where all three pairs have the invariant mass of 0(Mg), is power-suppressed. 
This contribution is the one, which can be calculated perturbatively in terms of the diagrams with two hard gluon 
exchanges. 

We define the power counting rules for the various topologies in the dominant kinematic region, and identify the 
leading ones. Consider first nonresonant contributions. Topology I behaves like (AMb)~^, where one power of 
(AMb)^^ comes from the hard gluon in Fig. 1(a), which kicks the soft spectator in the B meson into a fast one in a 
light meson [7], and another power is attributed to the invariant mass of the light meson pair. The overall product 
of the meson decay constants is not shown explicitly. Topology II exhibits the same power behavior as topology I: 
the hard gluon in Fig. 1(b), i.e., the B meson transition form factor, gives a power of {AMb)~^, and the light-light 
form factor gives another power. The scalar vertex introduces an extra power uiq/Mb, mo being the chiral symmetry 
breaking scale, to topology lis. Topology III must involve large energy release for producing at least a pair of fast 
mesons with the invariant mass of 0(Mg). That is, it behaves like (AMb)~^M^^. Hence, we have the relative 
importance of the decay amplitudes, 

M:A^n:Mis:A^m = l:l:^:4, (2) 

indicating that topology III is negligible. For resonant contributions, we replace the power of (AAf^)^^ associated 
with the light meson pair by A~^. Therefore, Eq. (2) still holds. 

Take topology I for the K^-k^'k~ mode as an example, in which the B meson transit into a pair of pions. The 

7r+ and tt" mesons carry the momenta Pi and P2, respectively. The B meson momentum Pb, the total momentum 
of the two pions, P = Pi -|- P2 , and the kaon momentum P3 are chosen as 

Mr Mr Mr 

Pb = -^(1,1,0t), P=-^(l,r?,OT), P3 = -^(0, 1 - ry, Or) , (3) 

with the variable ry = w'^ /Mg, = being the invariant mass of the two-pion system. The light-cone coordinates 
have been adopted here. Define ( = P^ /P^ as the 7r+ meson momentum fraction, in terms of which, the other 
kinematic variables are expressed as 

p^ = {i-Op^, pf = (1 - c)r?p+ , P2- = cr/p+, piT = piT = a^-(:W ■ (4) 

The two pions from the B meson transition possess the invariant mass uP' ^ O(AMb), implying the orders of 
magnitude P+ ~ 0{Mb), P~ ~ 0(A) and Pt ~ O(VAMb). In the heavy-quark limit, the hierachy P+ » Pi(2)t > 
P~ corresponds to a coUinear configuration. Therefore, we introduce the two-pion distribution amplitudes [13] , 

2^ e— ^^^"(7r+(Pi)7r-(P2)|V;(y-) /i-TV(0)|0) , (5) 

^^^"(7r+(Pi)7r-(P2)|Vi(y-)T^(0)|0) , (6) 

y~{7r+{P^)n-{P2my-)i(T^.n''_P-'Ti;{0)\0) , (7) 



with <I>t, being the twist-2 component, and $s and <I>t the twist-3 components. T = is for the isovector 1=1 state, 
tp the u-d doublet, z the momentum fraction carried by the spectator u quark, and n_ = (0, 1,0t) a dimensionless 
vector. The constant 72';^ of dimension of mass is defined via the local matrix element [13], 

hm (7r+(Pi)7r-(P2)|V5(0)2a^.n^P^^V'(0)|0) = |^(2C - 1)P+ • (8) 

w ^ J2-K 

The matrix element with the structure 75 vanishes for topologies I and lis, and contributes to topology II 
at twist 4. The one with the structure 75 vanishes. For topologies II and lis, a kaon-pion distribution amplitude 
is introduced in a similar way. For other two-pion systems, the distribution amplitudes can be defined with the 
appropriate choice of the matrix T. For instance, T = 1/2 is for the ir^n^ isoscalar (/ = 0) state. 
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A two-pion distribution amplitude can be related to the pion distribution amplitude through the calculation of the 
process 77* tt+tt" at large invariant mass [14]. The extraction of the two-pion distribution amplitudes from the 
B ttttW decay has been discussed in [15]. Here we pick up the leading term in the complete Gegenbauer expansion 
oi^i{z,C,w^) [13]: 



2, 3F^,t(u;2) ^^^^ ^ , ^ 2^ 3F,(u;2) 



-z(l-z)(2C-l), <^>s{z,C,w') 



/2n:. 



z{l-z) 



(9) 



where F7r,s,t(u>^) are the time-like pion electromagnetic, scalar and tensor form factors with F7r,s,i(0) = 1. That is, 
the two-pion distribution amplitudes are normalized to the time-like form factors. For $t in Eq. (9), we have adopted 
the parametrization, 



(7r+(Pi)7r-(P2)|V;(0)»a^.n^P''^V(0)|0) = ^ F,{w') {2C - 1)P+ . 

J 27r 

Note that the asymptotic functional form for the z dependence of is an assumption. 
For the B meson distribution amplitude, we employ the model [5], 



(10) 



$B(a;) = Nbx^{1 - xf exp 



1 fxMB\ 

'2 \ UB ) 



(11) 



with the shape parameter ub = 0.4 GeV, and the normalization constant Nb related to the decay constant /s = 190 
MeV (in the convention = 130 MeV) via /^^ <^B{x)dx — fB/{2y/2N^). The above $_b is identified as in the 
definition of the two leading-twist B meson distribution amplitudes $± given in [16,17]. Equation (11), vanishing at 
X ^ 0, is consistent with the behavior required by equations of motion [18]. Another distribution amplitude $s in 
our definition, identified as i's = {^b — ^b)/v^ with a zero normalization, contributes at the next-to- leading power 
A/Mb [7]. It has been verified numerically [19] that the contribution to the B ^ n form factor from $b is much 
larger than from $ b ■ 

The total decay rate is written as 



5127r4 



/ rf??(l -v) [ 
Jo Jo 



with the amplitudes. 



Ml 



v: 



i=4,6 



Mi = (2C-l)i^.(a;') 



v: 
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M = Mi + Mil + A^iis 



(12) 



Mn. = V*F,{uj')F:e 



p{d) 



i=3 



i=3,5 



(13) 



For a simpler presentation, we have assumed that the kaon-pion time-like form factor in topology lis is equal to 
the pion time-like form factor multiplied by the ratio of the decay constants ]kI This assumption is in fact not 
necessary, and the property of the kaon-pion form factor will be discussed elsewhere. The siiperscript P(g) stands for 
the amplitude from a penguin operator producing a pair of quarks q. Those without P(q') arise from tree operators. 
The subscript ei stands for the emission topology (in contrast to the annihilation topology III) from the effective 
four-fermion operator Oi in the standard notation. 

We calculate the hard kernels by contracting the structures, which follow Eqs. (5)-(7), 



{/Pb + Mb)75 



^b{x) , 



1 



w(2C - 1) 



V2Nc 

to Fig. 1. The factorization formulas for the B ^ 2tt transition amplitudes are given by 

^b{xi) 



(14) 



pi 

jr^^(") = SttCfMKi - ?7) / dxidz 

Jo 



xizMl + P} 



:{[{! + z)^y{z, C,w^) + ^{1 - 2z)Mz, C, «;') + y/v{l - 2^)$«(^, C, w"") 



zM% + P2 



3 



rj^y {z, C, w^) - 2 {z, C, w^) 



xiMl 



pi 

^^^(") = -IQ-KCpMlro I dxidz 
Jo 

:{[0 



+ P2 



(15) 



with To = mo/ Mb- Te.i is the same as T^^^^ but with 04" •* replaced by (here a2 is close to unity). The definitions 

of the Wilson coefficients a^'i^[t) arc referred to [20]. The hard scales are defined by t^i^^ = max[\/zAfB, Py] and 

t^^ = in.a.x[,yxi M B, Pt]- The above coUinear factorization formulas are well-defined, since the invariant mass of the 
two-pion system, proportional to Pt, smears the end-point singularities from z ^ 0. 
The B meson transition form factors involved in topologies II and lis are 

pl poo 

F^}"^^ = ^-kCfMI I dxidxs bidbib3db3<^B{xi,bi) 
Jo Jo 

X { [(1 - 77)(1 + (1 - v)x3)^k{x3) + ro(l + ry - 2(1 - ?7)a;3)*^(x3) 
+ro(l - 77)(1 - 2x3Wk{x3)] Ef\t^p)h,{xu (1 - 77)2:3, 61, 63) 
+2ro(l - nWK{xi)E'f\tf'^)he{{l - r])x3,xub3, b,)} , 

pl poo 

l^TTCpMB^/rj I dxidxs / 61^6163^63 $B(a^i,6i) 
Jo Jo 

X { [(1 - v)^k{x3) + 2ro$^(a;3) + ro(l - 77)a;3 {^^M - ^'ki^s)) 

xEf\t^^))he{x^,{l-7l)X3M,h3) 

+2ro(l - 77)$^(x3)4''^(4'^)/ie((l - ri)x3, a;i, 63, 61)} . 



^e6 



(16) 



The definitions of the evolution factors E\''\t), which contain the Wilson coefiicients a}f\t), of the hard functions 
he, and of the kaon distribution amplitudes <^k, ^^'^ referred to [20]. Pg^'^\ -^e^^*^^ ^^^'^ -^ei are obtained 

from F^^^'^^ by substituting 03'^ a^^ and ai for 04^*^ respectively. 

The PQCD evaluation of the form factors indicates the power behavior in the asymptotic region, Ft^Ivo^) ^ l/w)^, 
and their relative importance: Fg^tiw"^) / F-Kiw"^) ~ mo/w. Therefore, the twist-3 contributions in Eq. (15) are down 
by a power of ,Jrjmo/w = mo /Mb compared to the twist-2 ones, which is the accuracy considered here. To calculate 
the nonresonant contribution, we propose the parametrization for the whole ragne of w^. 



(9)/ 



Fi-^\w') 



m 



Es,t (w ) 



yj2 _|_ ^2 



mom 



+ mom^ ' 



(17) 



where the the parameter m = 1 GeV is determined by the fit to the experimental data M'j^_^\Fjr{M'j^^)\'^ ^ 0.9 GeV^ 
[21], Mjj^ being the J /ip meson mass. These form factors can carry strong phases, which are assumed to be not very 
different, i.e., overall and negligible here. 

To calculate the resonant contribution, we parametrize it into the time-like form factors, 



Ml 



Ml 



V(w2 - MIY + r^u;2 u;2 + Ml 



2 ' 



(18) 



with Fy being the width of the meson V . The subtraction term renders Eq. (18) exhibit the features of resonant 
contributions: the normalization Pi'^^(O) = and the asymptotic behavior Fi^\vP) ~ 1/w^, which decreases at 
large w faster than the nonresonant parametrization in Eq. (17). Equation (18) is motivated by the pion time-like 
form factor measured at the p resonance [22] . It is likely that all P^'^ t contain the similar resonant contributions. 
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The relative phases among different resonances will be discussed elsewhere by employing the more sophisticated 
parametrization [23]. Here we assume the absence of the interference effect. 

We adopt mo = 1.4 (1.7) GcV for the pion (kaon) and the unitarity angle 03 = 90° [5]. For the B+ p°{770)K+ 
and /o(980)iir+ channels, we choose Tp = 150 MeV and T/q =50 MeV [24]. The nonresonant contribution 

0.61 X 10~^ to the i^'+TT+TT" branching ratio is obtained. Our results 1.8 x 10~^ and 13.2 x 10~® are consistent 

with the measured three-body decay branching ratios through the — !■ p{770)K~^ and S+ — > /o(980)-fir"'" channels, 
< 12 x 10~^ and (9.6l2!3-i'.5-o.8) ^ 10~^ [11], respectively. Since the /o width has a large uncertainty, we also consider 
F/q =60 MeV, and the branching ratio reduces to 10.5 x 10~^. The resonant contributions from the other channels 
can be analyzed in a similar way. For example, the Jsr*(892) resonance can be included into the K-ir form factors 
by choosing the width F^* = 50 MeV. The nonresonant and resonant contributions to the — » if+TT+Tr" decay 
spectrum are displayed in Fig. 2. 

In the above formalism nonfactorizable contributions arise from the diagrams, in which a hard ghion attaches 
the spectator quark and the meson emitted from the weak vertex (topology I) or the meson pair (topologies II and 
lis). The nonfactorizable contributions, suppressed by ln~^(MB/A) [8], and topology III, being of 0(A/Mb), can 
be evaluated systematically by means of the two-meson distribution amplitudes. The framework presented here is 
not only applicable to the study of three-body mesonic B meson decays, but also to baryonic decays [25], such as 
B ppK. One simply introduces two-proton distribution amplitudes, and the calculation of the corresponding hard 
kernel is similar. 

In this letter we have proposed a promising formalism for three-body nonleptonic B meson decays. This formalism, 
though at its early stage, is general enough for evaluating both nonresonant and resonant contributions to various 
modes, and as simple as that for two-body decays. In the future we shall discuss more delicate issues, such as CP 
asymmetries [26], phase shifts from meson-meson scattering [27], and interference effects among different resonances 
[28]. 
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part by the National Science Council of R.O.C. under Grant No. NSC-91-2112-M-001-053, by the National Center 
for Theoretical Sciences of R.O.C, and by Theory Group of KEK, Japan. 
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FIG. 1. Graphic definitions for topologies I, II(s), and IIL 
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FIG. 2. (a) [(b)] Nonresonant (resonant) contribution to the K^-k^-k~ decay spectrum with respect to the two-pion 

invariant mass M(7r^7r~). The sharp peak corresponds to the /o resonance with the width 50 MeV. 
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